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Abstract 


We extend the parameter space of minimal supergravity to negative values of mg, the univer¬ 
sal scalar mass parameter defined at the grand unified scale. After evolving to the weak scale, 
all scalars can be non-tachyonic with masses consistent with collider constraints. This region of 
parameter space is typically considered excluded by searches for charged dark matter, since the 
lightest standard model superpartner is a charged slepton. However, if the gravitino is the lightest 
supersymmetric particle, the charged slepton decays, and this region is allowed. This region pro¬ 
vides qualitatively new possibilities for minimal supergravity, including spectra with light sleptons 
and very heavy squarks, and models in which the lightest slepton is the selectron. We show that 
the mg < 0 region is consistent with low energy precision data and discuss its implications for 
particle colliders. These models may provide signals of supersymmetry in even the first year of 
operation at the Large Hadron Collider. 

PACS numbers: 04.65.-|-e, 12.60.Jv, 13.85.-t 


1 


I. INTRODUCTION 


Supersymmetric models are theoretically motivated extensions of the standard model 
(SM) of particle physics that predict both direct and indirect signals in particle physics 
experiments. Most analyses of supersymmetric models assume the minimal supersymmet¬ 
ric standard model (MSSM), the supersymmetric extension of the SM that contains the 
minimal number of superpartners. Supersymmetry (SUSY) must also be broken. To make 
phenomenological analyses tractable, a moderately simple model for soft supersymmetry¬ 
breaking terms must be chosen. 

By far the most studied supersymmetric model is minimal supergravity (mSUGRA) [1], 
which is specihed by 6 parameters: 

ml, Mi/ 2 , Aq, tan(3, sign(/i), and mg/a . (1) 

Here ml is the universal soft scalar mass squared, Mi/a is the universal soft gaugino mass, Aq 
is the universal soft trilinear term, tan (3 is the ratio of the vacuum expectation values of the 
up and down type Higgs bosons, /i is the supersymmetric Higgs mass parameter, and ms/a 
is the gravitino mass. The hrst three terms are dehned at the grand unihed theory (GUT) 
scale Mgut — 2.4 x 10^® GeV, where the gauge couplings unify. All superpartner masses 
and couplings are determined by these 6 parameters and renormalization group equations 
(RGBs). The lightest SM superpartner is typically either the lighter stau or the lightest 
neutralino. 

Note that mSUGRA is typically thought to be determined by the first 5 parameters. 
When the gravitino is not the lightest supersymmetric particle (LSP), much of cosmology 
and all of particle phenomenology is insensitive to 7713/2 • However, if the gravitino is the 
LSP, both cosmology and particle phenomenology are sensitive to the gravitino mass, and 
7773/2 is an essential parameter of mSUGRA. 

When i?-parity is conserved, as we assume throughout this study, the LSP is stable. 
Gommonly it is (implicitly) assumed that the gravitino is not the LSP. In this case, the 
region of parameter space in which the stau is the lightest SM superpartner is strongly 
disfavored, as it predicts an absolutely stable charged massive particle (GHAMP), which has 
not been fonnd [2, 3]. Results of mSUGRA analyses are often displayed in the {ml,Mi/ 2 ) 
plane. Null results from GHAMP searches then exclnde from consideration a thin triangnlar 
wedge in this plane with small ml > 0 and, a fortiori, the entire half plane with ml < 0. 

As emphasized above, however, this line of reasoning relies on the assnmption that the 
gravitino is not the LSP. There are no theoretical motivations for this assumption — the 
gravitino mass is a free parameter in mSUGRA. In this, as well as in other scenarios with 
high-scale snpersymmetry breaking, it is naturally of the same order of magnitude as other 
superpartner masses, but it is not necessarily larger. In addition, recent work has estab¬ 
lished that there are also no phenomenological or cosmological reasons to exclude the grav¬ 
itino LSP scenario [4-11]. In fact, high-scale snpersymmetry breaking with a gravitino 
LSP has a number of novel implications and virtnes. For example, if the gravitino is the 
LSP and the next-lightest supersymmetric particle (NLSP) is charged, the signal of SUSY 
at colliders will be metastable charged particles. Such particles have spectacular signa¬ 
tures [12-15]. They also make possible the investigation of gravitational interactions and 
the qnantitative verification of snpergravity in high energy physics experiments [16-20]. Gos- 
mologically, gravitinos produced through decays naturally have the correct relic density to 
be snperweakly-interacting massive particle (snperWIMP) dark matter. For some parame¬ 
ters, gravitino dark matter produced in this way has featnres usnally associated with warm 
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dark matter and may resolve controversial discrepancies in halo profiles and the formation of 
small scale structure [21-25]. Last, the late decays that produce gravitino dark matter also 
release electromagnetic and hadronic energy, with (possibly felicitous) implications for Big 
Bang nucleosynthesis [4, 6, 9, 11, 26-29] and the cosmic microwave background [4, 6, 30]. 

In this paper we consider the possibility of mSUGRA with a gravitino LSP and ttIq < 0. 
We dehne 

mo = sign{ml)^J\ml\ . ( 2 ) 

We show that the region with mo < 0 contains models consistent with all collider limits. 
We also consider precision measurements, analyzing the anomalous magnetic moment of 
the muon a^, B(h ^ sj), and We hnd that the current discrepancy in 

between experiment and the SM prediction may be resolved for mo < 0 without disrupting 
the agreement for B{b ^ sj), and near future probes of B^ will have signihcant 

reach in mo < 0 model space. Precision data do not currently favor one sign of mo over the 
other. 

The simple modihcation of taking mo < 0 therefore “doubles” the viable mSUGRA 
parameter space and leads to qualitatively new possibilities. For example, in some regions 
of parameter space, the NLSP is not the stau, but the selectron. This overturns the common 
lore that Yukawa couplings in RGBs lower soft masses; when some scalars are tachyonic in 
part of the RG evolution, Yukawa terms may increase scalar masses. We also hnd that light 
charged sleptons can be produced for any value of Mi/ 2 - This produces spectra where the 
charged sleptons have masses around 100 GeV, but all other superpartner masses are above 
1 TeV, with squark and gluino masses around 3 — 4 TeV. Such spectra are not found for 
mo > 0 and have novel implications for the Large Hadron Gollider (LHG) and International 
Linear Gollider (ILG). Although all of these features may be found in general MSSM models, 
it is striking that we hnd them here in a framework with universal scalar and gaugino masses, 
motivated as these features are by simplicity, the SUSY havor and GP problems, and gauge 
coupling unihcation. 

Gosmologically, these models diher from conventional models with mo > 0 in several 
ways. As noted above, the mo < 0 models have superWIMP dark matter, as opposed to 
the conventional neutralino WIMP dark matter, with the implications mentioned above for 
Big Bang nucleosynthesis, the cosmic microwave background, and structure formation. In 
addition, there are possibly novel implications for vacuum stability [11] and gauge symmetry 
breaking at high temperatures. We defer discussion of cosmological issues [31], and focus 
here on implications for particle physics. 

This study is organized as follows. In Sec. II we show that, even given mo < 0 at the GUT 
scale, all superpartner masses, when evolved to the weak scale, can have values consistent 
with current experimental bounds. We delineate the allowed regions and determine which 
regions of parameter space have stau and selectron NLSPs. The resulting superpartner 
masses in the mo < 0 region are discussed in Sec. III. Low-energy observables are analyzed 
in Sec. IV. In Sec. V, we show two representative superpartner spectra and briefly discuss 
the implications for the LHG and ILG. In Sec. VI, we conclude and indicate interesting 
avenues for further investigation. 
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FIG. 1: Regions of the (mo,Mi/ 2 ) plane, extended to mo < 0, for Aq = 0, /r > 0, and tan/3 = 10 
(left) and tan/3 = 60 (right). The green (medinm shaded) region is experimentally excluded, and 
the unshaded region is the conventional neutralino (N)LSP region. In the remaining regions, the 
gravitino is the LSP: in the yellow (light shaded) region, the stau is the NLSP, and in the thin 
magenta (dark shaded) region of the tan /3 = 60 plot, the selectron is the NLSP. The present Higgs 
mass bound rrih > 114.1 GeV excludes regions below the solid contours. The symbols © mark the 
location of benchmark Models A and B; their RG evolution is shown in Fig. 2. 

II. REGIONS OF MSUGRA PARAMETER SPACE FOR mo < 0 

In this section, we determine the allowed regions of mSUGRA parameter space with 
mo < 0 , and further classify the allowed parameter space according to what particles are 
the LSP and NLSP, since these play a large role in determining experimental signatures. 

For mo < 0, an immediate worry is that scalar masses will remain tachyonic even after RG 
evolution to the weak scale. As usual, gauge interactions raise the soft masses, and so the 
most problematic scalars are the right-handed sleptons, since these have only hypercharge 
interactions. The RGBs for conventional mSUGRA have been studied in great detail. A 
well-known approximate relation for the weak-scale right-handed selectron mass in terms of 
GUT-scale parameters is [37] 

+ . (3) 

This remains valid for mo < 0. From this, we see that negative m^ can always be compen¬ 
sated by large M 1/2 to make the right-handed sleptons, and with them the entire superpart¬ 
ner spectrum, non-tachyonic. 

Allowed regions of the (mo,Mi/ 2 ) plane are shown for two values of tan/3 in Fig. 1. 
The SUSY spectra have been calculated with the software package ISAJET v7.71 [32], 
modihed to accommodate mo < 0. ISAJET includes 2-loop RGBs and 1-loop corrections 
to superpartner masses, and we choose a top quark mass of 175 GeV. The green (medium 
shaded) region is excluded. For mo > 0 the boundary is determined by the LEP chargino 
mass limit rny,± > 103.5 GeV [33]. For mo < 0, it is essentially determined by null searches 
for long-lived charged sleptons at LEP, leading to limits mj^ > 99 GeV [34]. For tan/3 = 10, 
the border for mo < 0 follows to a reasonable approximation the tachyonic slepton line 
mo = — 2 . 6 M 1/2 one can derive from Eq. (3). For tan/3 = 60, the excluded region has a 
more complicated shape because large 1 -loop corrections play an important role, as discussed 
below. 
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The allowed regions may be further divided according to what particles are the LSP and 
NLSP. The unshaded regions of Fig. 1 are the conventional regions in which the lightest SM 
superpartner is the lightest neutralino. It is either the LSP or, if the gravitino is the LSP, 
the NLSP. 

In the rest of the allowed regions shown, the gravitino must be the LSP to avoid having 
charged dark matter, and the NLSP is a charged slepton. To determine which charged 
slepton is the NLSP, consider the RGBs for their soft mass parameters. At 1-loop, these are 



where t = hi . As is well-known, when all mass parameters are non-tachyonic, 

Yukawa interactions lower soft masses, leading to the common lore that, given universal 
scalar boundary conditions, the lightest slepton is always the stau. In the present case, 
however, mo < 0, and so in evolving from the GUT scale, selectron masses initially rise 
slower than stau masses. Of course, for the spectrum to be viable, all physical scalar 
masses must eventually become non-tachyonic, and so will exert the conventional effect of 
Yukawa couplings as one approaches the weak scale. (The Higgs scalar mass parameters may 
remain negative.) The competition between the new and the conventional effects determines 
whether the NLSP is the selectron or the stau. 

These effects are shown in Fig. 2 for the two benchmark models highlighted in Fig. 1. 
In the left panel, all scalar masses of Model A quickly become positive as they evolve from 
the GUT scale, and so the stau becomes the NLSP, as usual. In the right panel, however, 
the scalar masses of Model B are negative for much of the RG evolution, and becomes 
negative, leading to an inverted flavor spectrum with a selectron NLSP. 

In Fig. 1, the fi is the NLSP in the yellow (light shaded) region, and cr is the NLSP in 
the magenta (dark shaded) region. The current experimental limits force sleptons to be not 
just non-tachyonic, but signihcantly so, and so the viable selectron NLSP region is reduced 
to a thin sliver near the upper, left-hand corner in the tan/? = 60 plot. Its exact location 
is therefore rather sensitive to small corrections and depends on the implementation of RG 
evolution and loop-level corrections to the superpartner mass spectrum. The mere possibility 
that universal slepton masses can lead to selectrons lighter than staus, however, is a robust 
physics effect; it is novel and never realized in conventional mSUGRA. The selectron NLSP 
region may be much larger in even slightly more general models. For example, motivated 
by SO(IO) unihcation, one may consider models in which the matter scalar masses are 
determined by the parameter ttiq, but the Higgs scalar masses are unified at a different 
GUT-scale parameter rn\j. In these non-universal Higgs mass models [35], by choosing 
< ttIq < 0, the stau mass will receive large positive contributions from that are 
absent for the selectron, as can be seen in Eq. (5), and the selectron NLSP region will be 
much larger. Last, we note that in the mSUGRA models considered here, the possibility 
of a non-tachyonic sneutrino (N)LSP also exists at low M 1 / 2 , but this lies entirely in the 
excluded region. 

As noted above, the excluded region can have a rather complicated shape. For tan (3 = 60, 
for example, as can be seen in the right-hand panel of Fig. 1, the excluded region has an 
interesting shape. This results from the remarkable fact that the light charged slepton 
masses do not increase monotonically as one increases M 1/2 for constant tuq. In fact, this is 
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FIG. 2: The RG evolution of soft scalar masses in Model A, the stau NLSP point with uiq = 
—40 GeV, Mi /2 = 300 GeV, tan /? = 10 (left) and in Model B, the selectron NLSP point with 
mo = —700 GeV, M 1/2 = 1900 GeV, tan/3 = 60 (right). In both cases, Aq = 0 and /i > 0. 



FIG. 3: Regions of the (mo,Mi/ 2 ) plane, extended to mo < 0, for Ao = 0, /r > 0, and tan/3 = 60 
as in the right-hand panel of Fig. 1, but with 1-loop corrections to sparticle masses neglected. 


not peculiar to mo < 0, as it occurs even for positive constant mo. 

This behavior results from 1-loop corrections present in the slepton mass matrix 

( ^LL + ^LR + \ /cN 

where are tree-level contributions and 6M^ are 1-loop corrections. If 1-loop corrections 
are neglected, the slepton mass monotonically increases for increasing mo and fixed M 1/2 (or 
vice versa), as illustrated in Fig. 3. However, for large tan/3 the loop corrections have the 
proper sign and magnitude to lower the mass eigenvalues of the stau below the experimental 
bounds in part of the parameter space. Of course, the 1-loop corrections are physical, and 
we include them in all plots and results below. 

To conclude this section, the impression that neutralinos are the lightest SM superpart¬ 
ners in most of mSUGRA parameter space is artihcial: it follows only if one requires mo > 0. 
Allowing mo < 0 extends the viable region of mSUGRA parameter space signihcantly and 
shows that staus are the lightest SM superpartners in much of mSUGRA parameter space. 
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FIG. 4: Squark masses for mSUGRA extended to uiq < 0 for Aq = 0, /r > 0, and tan/3 = 10 (left) 
and tan/3 = 60 (right). The up-type squarks (top) are ul (solid black), uji (dotted black), ii (dot 
dashed black), and t 2 (dashed black). Similarly the down-type squarks (bottom) are (solid 
black), dpi (dotted black), bi (dot dashed black), and 62 (dashed black). The gluino mass (cyan, 
solid light) is presented on all 4 plots. The contours are for masses 1 TeV, 2 TeV, and 3 TeV from 
bottom to top. 


In addition, allowing mo < 0 leads to other new phenomena, such as the possibility that the 
selection is the lightest SM superpartner. In the next section, we explore the implications 
of mo < 0 for the sparticle spectrum more fully. 


III. SUSY MASS SPECTRA FOR mo < 0 

The squark and gluino masses are presented in the (mo, M 1 / 2 ) plane in Fig. 4. In the 
allowed region with mo < 0, they are relatively insensitive to mo, since they are dominated 
by the RG contributions of the gaugino masses. The left-handed down and up squarks are 
approximately degenerate, as are the right-handed up and down squarks. The gluino is 
always heavier than all squarks in the mo < 0 allowed region. 

Slepton masses are shown in Fig. 5. In contrast to squark masses, slepton masses are 
extremely sensitive to mo in the allowed tuq < 0 region. Contours of constant slepton mass 
switch from concave down for mo > 0 to concave up for mo <0. As a result, light sleptons 
near their experimental limit may be found for any value of M 1 / 2 . This is true even though 
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FIG. 5: Slepton masses for mSUGRA extended to uiq < 0 for Aq = 0, /r > 0, and tan/3 = 10 (left) 
and tan/3 = 60 (right). The contours are for fi (dashed), ei (solid) and cr (dotted). The contours 
are for masses 200 GeV, 500 GeV, and 1 TeV from bottom left to top right, with the exception 
that in the right panel, because the slepton masses drop and then rise again for fixed M 1/2 and 
increasing niQ, there are two sets of 200 GeV contours for both cr and fi. For both of these 
particles, the leftmost 200 GeV contour is barely visible in the selectron NLSP region. Throughout 
the parameter space, the sneutrinos and f '2 are almost degenerate with the cr. 



mo mo 

FIG. 6; Neutralino and chargino masses for mSUGRA extended to mg < 0 for Ag = 0, /r > 0, 
and tan/3 = 10 (left) and tan/3 = 60 (right). The lightest neutralino Xi (dashed) has values from 
bottom to top of 200 GeV, 500 GeV, and 800 GeV. The lightest chargino xf (dotted) has values 
from bottom to top of 500 GeV, 1 TeV, and 1.5 TeV, while /r (solid) has values 500 GeV, 1 TeV, 
and 2 TeV from bottom to top. 

the masses of all other superpartners becomes large for large M 1 / 2 , as may be seen in Model 
B, the selectron NLSP model shown in the right panel of Fig. 2. 

In Fig. 6, we present contours for neutralino and chargino masses and for the Higgsino 
mass parameter ft. In the mg < 0 region, \fi\ is always much larger than the electroweak 
gaugino masses Mi and M 2 . As a result, the lighter chargino and lighter two neutralinos 
are nearly pure gauginos, with Xi ~ -B, X2 ~ xt ~ W^, and m^o Mi and m^o 
m^± ~ M 2 ~ 2Mi. 

Last, Higgs boson masses are given in Fig. 7. The masses of the Higgs bosons, h^, A^, 
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FIG. 7: Higgs boson masses for mSUGRA extended to mo < 0 for Hq = 0, /i > 0, and tan/? = 
10 (left) and tan/3 = 60 (right). The SM-like Higgs boson hP mass contours (solid) have values 
114 GeV, 118 GeV, and 122 GeV from bottom to top. The AS' mass contours (dashed) have values 
500 GeV, 1 TeV, and 1.5 TeV from bottom to top. The heavy GP-even and charged Higgs scalars 
are approximately degenerate with the A^. 

and are all insensitive to ttiq for mo < 0. The predicted valne of ruho is above 114 GeV 
throngliont the allowed mo < 0 region, and so consistent with cnrrent bonnds. It increases 
with increasing M 1 / 2 , rising to approximately 122 GeV for M 1/2 = 1 TeV. and 

are all approximately degenerate for mo < 0. 




IV. PRECISION EXPERIMENTAL CONSTRAINTS 

We now consider constraints on these models from precision experimental data. We focns 
on three processes that are well-known to have signihcant sensitivity to snpersymmetric 
contribntions: the anomalons magnetic moment of the mnon and the rare decays b —>■ sj 
and These contribntions have been calcnlated nsing the software package 

micrOMEGAs, vl.3.6 [38]. 


A. Anomalous Magnetic Moment of the Muon 

Determining the SM valne of the anomalons magnetic moment of the mnon is not straight¬ 
forward becanse of hadronic contribntions to higher order loop processes. These hadronic 
loop contribntions are nsnally estimated from measnrements of e’''e“ —> hadrons or r —> 
hadrons. The resnlting SM predictions for are = 116 592 018 (63) x 10“^^ if the r data 
are nsed, and = 116 591835 (69) x 10“^^ if the data are nsed [39, 40]. Given theo¬ 
retical assnmptions reqnired to nse the r data, the valne is generally jndged to be more 

reliable [41]. These shonld be compared to the measnred valne = 116 592 080 (60) x 10“^^ 
from the Mnon {g — 2) Gollaboration [42, 43]. Taking the SM valne for nsing the 
data, there is a discrepancy between theory and experiment of = 245 x 10“^^, a deviation 
of approximately 3a. 

The SUSY contribntion to in mSUGRA with mo < 0 is shown in Fig. 8. A deviation 
consistent with the discrepancy between experiment and the SM prediction may be 
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FIG. 8; The SUSY contribution to for mSUGRA extended to tuq < 0 for Aq = 0, /r > 0, and 
tan/3 = 10 (left) and tan/3 = 60 (right). From bottom to top the contour values are 300, 250, 200, 
150, and 100 in units of 10“^^. 

obtained for /r > 0 and light nentralinos and sleptons or light charginos and snentrinos. As 
noted in Sec. Ill, sleptons are light along the entire mo < 0 experimentally exclnded border, 
bnt the ganginos increase in mass as M 1/2 increases. A large SUSY contribntion to is 
therefore fonnd only for relatively small M 1 / 2 . The 3a deviation mentioned above may be 
explained, for example, for tan/3 = 10 and M 1/2 ~ 300 GeV. 


B. B(b —^ S 7 ) 

The flavor changing nentral cnrrent transition b ^ sj has a branching fraction measnred 
to be 

r 3.21 ± 0.43 ± 0.27;[|-}g x 10“^ (CLEG) [44] 

B{b ^ 57 ) = i 3.88 ± 0.36 ± 0.3?;;];^^ x 10“^ (BABAR) [45] . (7) 

[ 3.55 ± 0.32 ± 0.30tE];J^ x 10“^ (BELLE) [46] 

The SM prediction is 3 . 79 to ;53 x 10“^ [47]. The mo < 0 mSUGRA predictions for B{b —>• 57 ) 
are shown in Fig. 9. The experimental and SM theory valnes agree within errors, and so 
b —>■ sj may, in principle, eliminate models with light sqnarks and ganginos. As can be 
seen in Fig. 9, however, snpersymmetric effects in the plotted mg < 0 region are never large 
enongh to create a discrepancy between these mSUGRA models and experiment. 

C. 

The branching fraction for B^ decaying to two leptons is an important measnrement for 
constraining snpersymmetric models with large tan/3 [48-52]. The decay is enhanced by 
(tan/3)® for large tan/3. The cnrrent experimental bonnd is B{B^ —>■ /i"*"//”) < 1.5 x 10“^ 
from GDF II, based on 364 pb“^ of data [53], while the SM prediction is 3.42(54) x 10“® [54]. 
For small tan/3, mSUGRA with mo < 0, along with other SUSY models, predicts deviations 
far below cnrrent experimental bonnds. These deviations will not be probed nntil the LHG. 
However, for large tan/3, observable deviations are predicted even in Tevatron data. As 
shown in Fig. 10, the cnrrent Tevatron data do not exclnde additional parameter space. 
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FIG. 9: B{b 57 ) for mSUGRA extended to mo < 0 for Aq = 0, ^ > 0, and tan/3 = 10 (left) and 
tan/3 = 60 (right). The values from bottom to top are 3.10, 3.25, 3.40, and 3.55 in units of 10“^. 
In the tan/3 = 10 panel, B{b 57 ) does not exceed 3.59 x 10“^. 



FIG. 10: B{Bg —f ) for mSUGRA extended to uiq < 0 for Aq = 0, /r > 0, and tan/3 = 60. 
The lower contour is GDF II’s experimental upper bound 1.5 x 10“^, and the upper contour is for 
R(R° ^ /U+Ai“) = 1.0 X 10-8. 


Nevertheless, future improvements to sensitivities of ~ lO”® will probe the mo < 0 region 
all the way up to M 1/2 ~ 1.8 TeV for tan/3 = 60, and will also be sensitive to models with 
more moderate values of tan /3. 


V. COLLIDER SIGNALS 

As is well-known, the LHC will provide an extremely powerful probe of weak-scale su¬ 
persymmetry in the next few years. Here we discuss the implications of mo < 0 mSUGRA 
for the LHC and the proposed ILC. 

The prototypical signature of (R-parity conserving) supersymmetry at hadron colliders 
is missing transverse energy. In the case of mSUGRA with mo < 0, however, all SUSY 
events result in the production of two metastable charged sleptons. These have lifetimes of 
seconds to months, and so pass through collider detectors without decaying. These models 
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FIG. 11: Superpartner spectra for Model A, the stau NLSP point mo = —40 GeV, M 1/2 = 300 GeV, 
tan/3 = 10 (left) and for Model B, the selectron NLSP point mo = —700 GeV, M 1/2 = 1900 GeV, 
tan/3 = 60 (right). In both cases, Aq = 0 and // > 0. These models correspond to the parameter 
points highlighted with 0 symbols in Fig. 1; the RGBs of their scalars are displayed in Fig. 2, and 
their weak-scale parameters are given in Table I. 


ISAJET Specification 

Parameter 

Model A 

Model B 

MSSMA 

rUg n 

m^o tan/3 

720.70 395.95 

436.09 10.00 

3964.95 1994.33 

1274.81 60.00 

MSSMB 


652.21 601.12 603.80 

199.04 103.81 

3352.61 3178.77 3198.46 

1027.84 67.57 

MSSMG 

"^93 

At Ab A^ 

578.76 598.21 502.38 

198.49 101.56 

548.08 805.29 186.50 

2987.96 2860.45 2727.32 

1040.40 255.73 

2788.73 3147.17 14.15 

MSSMD 


Same as MSSMB (default) 

Same as MSSMB (default) 

MSSME 

All M2 

120.17 231.47 

831.43 1527.66 


Sag 

296 X 10“^^ 

78.1 X 10“^^ 


B {b ^ 57 ) 

3.49 X 10"^ 

3.57 X 10“^ 


B (RO ^ 

3.21 X 10"® 

8.84 X 10“^ 


TABLE 1: Mass parameters in GeV and the predicted values for precision observables for benchmark 
Models A and B. The masses in category MSSMA are physical masses; all other masses listed are 
soft SUSY-breaking parameters specified at the electroweak scale. 


therefore provide a conventional setting for what might otherwise be considered to be rather 
exotic signals, such as highly ionizing tracks and time-of-flight signatures [12-15]. Even a 
few events will provide unmistakable signals. 

As examples, let us consider the benchmark models indicated in Fig. 1. The RGBs for 
scalars in these models were shown in Fig. 2. In Fig. 11, we display the full superpartner 
spectrum for each of these models, and in Table I we list all mass parameters, which dehne 
these models at the electroweak scale. 

In Model A, the stau NLSP model of the right panel in Fig. 11, all superpartners have 
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masses under 1 TeV. This model is an excellent benchmark model. It explains the 3a 
deviation in a^, and preserves the agreement between the SM and experimental values of 
B{b ^ s'j). In addition, electroweak symmetry is broken radiatively and naturally, with 
fi ~ 400 GeV. 

This stau NLSP model has a total SUSY cross section at the LHC of (Tlhc( 14 TeV) = 
13.6 pb, as determined by ISAJET [32]. Even with an integrated luminosity of 100 pb“^, 
this implies over 1000 SUSY events, each with two metastable sleptons. In many of these, 
the sleptons will be slow enough to be seen as highly-ionizing tracks, providing a spectacular 
signal of new physics in even the first year of LHC operation. With more luminosity, large 
numbers of sleptons may be collected and their decays studied, making possible a variety of 
measurements with implications for cosmology, astrophysics and supergravity [4, 6, 9, 11, 16- 
28, 30]. 

Model B, the selectron NLSP model indicated in Fig. 1, is a complementary benchmark 
model. As seen in Fig. 11, the model has squarks and gluinos around 3 — 4 TeV, neutralinos, 
sneutrinos, non-SM type Higgs bosons, and left-handed sleptons around 1—2 TeV, a 210 GeV 
stau, 160 GeV selectron and smuon, and hnally a 124 GeV Higgs boson. The squarks and 
gluinos are too heavy to be produced with large cross section at the LHC. The biggest sources 
of SUSY particles at the LHC are direct Drell-Yan production of NLSP pairs, leading to 
2 metastable charged sleptons, and Drell-Yan production of the heavier sleptons, leading 
to even more unusual events with 2 taus, 2 muons/electrons, and 2 metastable charged 
sleptons. The total SUSY production cross section for this model is (Tlhc( 14 TeV) = 41 fb. 
This signal will be challenging to hnd in the first year of LHC running, but at the target 
luminosity of 100 fb“^/yr, the LHC will produce 4100 SUSY events per year. 

At the ILC, the total SUSY production cross sections are (Tilc( 500 GeV) = 1.35 pb for the 
stau NLSP Model A and (Tilc( 500 GeV) = 137 fb for the selectron NLSP Model B. Although 
squarks and gluinos are out of reach, the ILC will produce large numbers of sleptons. All 
the usual advantages of the ILC will allow detailed studies of the SUSY parameter space. 
In addition, the ability to produce sleptons at low velocities implies that they may be more 
easily trapped and studied than at the LHC. 


VI. CONCLUSIONS 

In this study, we have extended the well-known framework of mSUGRA to mo < 0. To 
our knowledge, this part of parameter space has not been considered previously, perhaps 
because it contains a charged slepton as the lightest SM superpartner. If the gravitino is 
the LSP, however, cosmological difficulties with CHAMPs are avoided, and this extended 
parameter space is allowed. We have noted that it is consistent with all limits from direct 
searches and constraints from low-energy precision measurements. In addition, we hnd 
models with qualitatively novel mSUGRA superpartner spectra, which predict spectacular 
signals with metastable charged sleptons at future colliders. 

Some of the cosmology of these models will be considered in a future work. It would 
be interesting to extend this work to more general models. For example, as argued above, 
we expect that the selectron may emerge as the NLSP generically in models with non- 
unihed Higgs masses, and there may well be other interesting phenomena. It would also 
be worthwhile to determine the reach of the LHC for various luminosities in the mo < 0 
parameter space; given how spectacular the signal of metastable charged particles will be, 
these models provide a welcome example in which supersymmetry may be discovered and 
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studied in even the first year of LHC running. 
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